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FAST  ION  CONDUCTION  IN  BORATE  GLASSES  B203-xLi 20-yLi 2S04 


1.1.  Structure  of  borate  glasses 
1.1.1  Introduction 

Fast  ion  conducting  glasses  are  generally  boron  oxide  or 
boron  sulfur  glasses  doped  with  mobile  ions  such  as  Li+,  for 
example.  Fast  ion  conducting  glasses  can  be  regarded  as  a 
solution  of  a  dopant  salt  such  as  Li2S04  in  a  matrix  composed  of 
the  glass  forming  materials  such  as  B203  and  its  modifier  Li20. 
Thus,  a  glass  such  as  B203-xLi20-yLi2S04  can  be  considered  as  a 
solution  where  Li2SC>4  is  the  solute  and  B203-Li20  is  the  solvent. 

Investigating  the  molecular  vibrations  of  the  different 
components  by  Raman  and  infrared  spectroscopies  gives  the 
possibility  of  performing  a  structural  analysis  and  eventually  of 
building  up  an  overall  representation  of  the  glass  structure. 

We  shall  examine  first  the  evidence  for  the  existence  of  the 
different  building  blocks  and  then  we  shall  concentrate  on  the 
possibility  of  inducing  structural  changes  by  varying  external 
physical  parameters  such  as  the  temperature,  for  example. 

1.1.2.  Structure  of  the  glass  former  B203 

The  borate  glasses  v-B203  are  generally  good  ionic  and 
electronic  insulators.  There  are  strong  arguments,  from 
spectroscopic  data,  that  the  structure  contains  mainly  rings 
whose  breathing  mode  is  Raman  active  and  shows  up  as  a  sharp  peak 
at  806  cm”l.  The  existence  of  such  a  sharp  structure  should  be 
surprising  for  a  glass  where  one  should  expect  to  observe  only 
very  broad  bands.  The  reason  is  that  the  sharp  structures  are 
due  to  highly  regular  hexagonal  rings  of  bonds  connected  into  a 
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more  disordered  network.  It  is  supposed  that  the  vibrational 
mode,  at  806  cm-1,  of  these  rings  is  decoupled  from  the  general 
network  and  consequently  is  strongly  localized.  A  condition  for 
vibrational  decoupling  of  the  appropriate  ring  mode  from  the 
surrounding  network  has  been  derived  ^ .  This  decoupling  is 
accomplished  by  a  near  cancellation  of  central  and  noncentral 
forces  within  the  ring  acting  on  the  high-coordination  B  atoms 
that  connect  to  the  network.  Therefore  the  connecting  atoms  do 
not  move  during  the  ring  vibration,  and  this  explains  the  absence 
of  a  mode  isotope  shift  when  isotopic  substitution  is  carried  out 
on  the  connecting  3  atoms.  A  representation  of  the  structural 
configuration  of  V-B2O3  is  a  given  in  Figure  1. 

1.1.3.  Effect  of  the  modifier  Li20 

The  addition  of  Li20  to  the  B2O3  modifies  the  structure  of 
the  glass.  As  the  oxygen  content  increases  in  the  glass  with 
increasing  concentration  of  Li2Or  an  increasing  number  of  B  atoms 
becomes  tetrahedrally  coordinated.  The  maximum  number  of  boron 
atoms  with  coordination  number  four  (BO^)  occurs  for  a 
concentration  of  x  of  Li20  of  the  order  of  40%^.  The  addition 
of  the  modifier,  Li20,  into  the  glass  B2O3  leads  to  the 
dissociation  of  the  metal  oxide  with  the  formation  of  a 
tetrahedral  radical  and  metallic  ions 
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A  BO4  unit  can  be  introduced  into  a  boroxol  ring  in  various 
ways (3) # 
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V~  B2°3 


Figure  1.  Structure  of  B2O3:  (a)  triangle  of  BO3;  (b)  planar 

ring  of  boroxol,  stabilized  by  delocalized  it  bonds; 

(c)  representation  of  a  random  network;  (d)  structure 
of  B2O3  with  rings  of  boroxol. 
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The  simplest  hypothesis  is  to  imagine  that  the  boroxol  rings  are 
progressively  replaced  by  triborate  groups  (Figure  2) . 

From  conductivity  measurements^1  one  can  deduce  that 
B203-xLi20  is  a  purely  ionic  conductor  with  Li+  as  the  mobile 
species . 

Increasing  the  L^O  concentration  x  causes  the  frequency  of 
the  breathing  mode  attributed  to  the  630^  ring  with  B04  groups  to 
decrease,  w  =  770  cm  ^ ,  and  the  asymmetry  of  the  band  to  increase 
toward  the  lower  energies  as  shown  in  Figure  3.  This  effect 
seems  to  be  moderated  by  the  addition  of  Li2S04. 

These  changes  of  frequencies  and  half-width  may  be  due  to 
the  appearance  in  the  matrix  conformation  of  more  complex 
configurations  in  the  8305  rings  such  as  tetraborate  and 
pentaborate  groups  shown  in  Figure  4. 

1.1.4.  Effect  of  the  dopant  L^SO^ 

The  addition  of  Li2S04  up  to  the  concentration  used  in  the 
experiments,  y  =  0.42,  does  not  seem  to  significantly  affect  the 
vibrational  frequencies  of  the  borate  matrix.  The  matrix  is 
formed  by  8305  rings  containing  one  or  more  B  atoms  of 
coordination  number  4,  whose  characteristic  frequency  is  780 
cm-1-.  In  the  case  of  borate  glasses  with  x  =  0.5  the  rings  are 
probably  mainly  triborates  but  the  existence  of  diborates  and 
ditriborates  cannot  be  excluded. 

The  presence  of  the  dopant  Li2S04  in  lithium  borate  glasses 
can  be  easily  identified  by  Raman  spectroscopy. 

The  Raman  spectrum  of  B2O3-O . 7Li 2^-0 . 42Li2S04  shown  in 
Figure  5  is  clearly  the  superposition  of  the  Raman  spectrum  of 
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Figure  2.  Some  of  the  structural  groups  appearing  in  the  alkal 
borates:  (a)  triborate?  b)  ditriborate?  c)  diborate 


RAMAN  INTENSITY  (orb.  units) 


700  780  800  850 

FREQUENCY  (cm-1) 


Figure  3.  Raman  peaks  corresponding  to  the  breathing  symmetrical 
mode  of  the  8305  cycles  as  a  function  of  the  Li 
concentration . 

Ai  B2O3-O . 5Li 2O  and  B2O3“0.7Li2O 

B:  B2O3-O . 5Li 2O-0 . 15Li 2SO4  and  B2O3-0.7Li2O-0.42Li2SO4 
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tetraborate  group 


pentaborate  group 


Figure  4.  Structural  group  appearing  in  the  alkali  borates: 


tetraborates  and  pentaborates . 


100  300  500  700  900  MOO 

RAMAN  SHIFT  (cm"') 


Figure  5 


Raman  spectra  for  B2O3-O . TL^O-y-L^SO^  for  y»0;  0.42 
at  room  temperature  without  annealing.  Laser  beam  at 
514.5  nm,  resolution  2  cm“^. 


Peaks 


the  matrix  B2O3-0.7Li2O  and  that  of  the  solute  Li2S04. 
such  as  v^,  and  are  those  of  the  characteristic 

vibrational  modes  of  S0~“ .  The  sharp  peak  shown  on  top  of  Figure 
1  is  that  of  the  matrix  B2C>3-0 . 7Li  20 . 

1.1.5.  Temperature  annealing,  glass  crystalline  transition  in 
borate  glasses 

Progressive  thermal  annealing ^  of  borate  glasses  leads  to 
the  crystalline  phase  of  B2°3-xL* 220-yLi 2S04  through  successive 
configurational  transformations  of  the  matrix.  The 
experiments^  are  conducted  as  follows:  the  samples  are 
submitted  by  slow  heating  to  a  given  annealing  temperature  Ta  and 
held  at  that  temperature  for  one  hour;  then,  are  rapidly  quenched 
to  room  temperature.  After  this  treatment,  the  Raman  spectra  are 
taken  at  room  temperature.  We  will  examine  successively  the 
different  parts  of  these  spectra. 

1.1. 5.1.  The  breathing  mode  v^(A^)  of  SO^ 

When  the  annealing  temperature  increases  progressively  from 
470°C  to  550°C  the  characteristic  peak  for  the  S0~“  mode  situated 
initially  at  1006  cm“^  first  decreases  in  intensity  and  then 
disappears.  A  new  peak  then  appears  at  1016  cirT^  which  has  still 
the  frequency  of  the  S04”  breathing  mode  but  now  in  the 
crystalline  phase  of  Li2S04.  This  transformation  can  be  followed 
closely  in  Figure  6. 

The  temperature  at  which  the  crystallization  of  Li2S04 
occurs,  signalled  by  the  appearance  of  the  peak  at  higher 
frequencies  near  1016  to  1020  cm-^/  depends  on  the  borate 
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RAMAN  INTENSITY  (orb.  units) 


Figure  6.  Raman  spectra  for  different  annealing  temperatures  of 
82O3-O. 5Li20-0. 15Li2S04  centered  at  the  mode  v^(A^). 
Spectra  taken  at  room  temperature  with  a  laser  beam  at 
514.5  nm  having  a  power  of  250  mW  and  with  an 
instrumental  resolution  of  2  cm"^. 
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composition.  For  B2O3-G . 5Li 2O-O . 1L i 2SO4  the  new  peak  at  1016 
cm"1  shows  up  at  516°C  and  is  fully  developed  after  the 
d isappea ranee  of  the  peak  at  1006  cm"1  at  540°C.  For 
^203-0 . 5L i 2O-O . 15Li 2S04  the  peak  at  1016  cm"^  appears  at  525°C 
and  the  1008  cm"1  peak  disappears  completely  at  579°C;  whereas 
for  B2O3-O . 7Li 2O-O . 42Li 2SO4  the  second  peak  is  already  clearly 
visible  at  475°C  and  shows  a  frequency  of  1020  cm-1.  The  first 
peak  initially  situated  at  1012  cm-1  disappears  above  511°C.  The 
remark  to  be  made  here  is  that  the  higher  the  concentration  of 
Li20  and  L^SO^,  the  lower  the  annealing  temperature  at  which  the 
Ll2f504  crystallites  appear  and  the  stiffer  the  SO^“  mode. 

The  stiffening  of  the  breathing  mode  probably  occurs  through 
electrostatic  interactions.  The  higher  the  ion  concentration  is, 
the  stronger  the  repulsive  force  is  among  the  SO^~  ions.  This  is 
an  effect  equivalent  to  increasing  pressure  and  therefore 
stiffening  the  symmetric  mode  frequencies. 

1.1. 5. 2.  The  breathing  mode  of  the  B3O5  ring 

In  the  glasses  with  the  diborate  matrix  B2O3-0. 21,^0  above 
an  annealing  temperature  of  510°C  one  observes  the  appearance  of 
a  new  peak  at  725  cm-1,  as  shown  in  Figure  7. 

For  annealing  temperatures  above  560°C  the  observed  spectrum 
is  identical  so  that  of  crystalline  diborate  B2O3-O . 5Li 2O ^ . 

The  crystallization  of  the  borate  matrix  into  the  forms  of 
diborate  is  also  accompanied  by  the  appearance  of  a  peak  at  1030 
cm-1  at  525°C  as  evidenced  by  spectra  shown  in  Figure  6. 

One  should  notice  also  the  appearance  in  the  crystalline 
phase  of  a  weak  peak  situated  at  772-774  cm"1  attributed  to  the 
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RAMAN  INTENSITY  (orb.  units) 


Figure  7.  Raman  spectra  of  B2O3-O . SL^O-O . lSL^SO^  for  different 
annealing  temperatures  focussed  at  the  frequency  region 
of  the  characteristic  modes  of  crystalline  diborate. 
Spectra  taken  with  a  laser  beam  of  514.5  nm  and  power  of 
250  mW.  The  instrumental  resolution  is  2  cm”1. 


/  Q  \ 

formation  of  some  ditriborate  groups' 

Depending  on  the  glass  matrix  composition,  the 
crystallization  of  the  matrix  may  occur  at  a  lower  temperature 
than  the  dopant  salt  if  x  <  0.5  and  after  the  crystallization  of 
the  salt  if  x  >  0.5.  The  respective  crystallization  temperatures 
for  different  glass  compositions  are  listed  in  Table  1.^)  T^ 
and  Tg  are  the  crystallization  temperatures  of  the  borate  glass 
matrix  and  dopant  salt,  respectively. 


TABLE  1 

CRYSTALLIZATION  TEMPERATURES  FOR  DOPED  BORATE  GLASSES 


B  20  3  —  x  L  i  20- 

-yLi 2SO4 

Tb  (°  C) 

Tg  (°  C) 

X 

y 

0.5 

0.05 

514 

529 

0.5 

0.10 

516 

536 

0.5 

0.15 

519 

525 

0.7 

0.2 

525 

478 

0.7 

0.4 

522 

500 

0.7 

0.6 

526 

518 
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1.2 


Ionic  conductivity  in  lithium  borate  glasses 


1.2.1.  Introduction 

The  electrical  properties  of  I^C^-xL^O-yL^SO^  have  been 
studied  by  the  method  of  complex  impedance  in  the  range  of 
10~3  -  10-4  Hz.  Gold  evaporated  electrodes  have  been  used  and 
the  measurements  performed  in  the  temperature  range  between  20 
and  400°C.  Determination  of  the  transport  number  by  the  method 
of  blocking  electrodes  ■ 9 )  confirms  that  the  material  is  an 
electronic  insulator. 

Figure  8  shows  the  variation  of  the  conductivity  as  a 
function  of  the  temperature  ^ 3-03  .  When  x  and  y  are  relatively 
high  the  conductivity  is  good,  o  =  2.10  Q  cm-  at  300°C,  and 
the  activation  energy  is  of  the  order  of  the  0.58  eV. 

Figure  9  gives  the  evolution  of  the  conductivity  at  300°C  as 
a  function  of  the  Li2SO^  concentration  y  for  the  different 
content  of  Li20,  x.  For  a  given  x  the  conductivity  increases 
with  increasing  y.  The  conductivity  rise  is  very  rapid  for 
x  <  0.43  and  becomes  much  slower  for  x  >  0.71.  A  saturation 
occurs  which  corresponds  to  the  non-dissociation  of  Li2SO^  and 
its  crystallization  which  limits  the  vitreous  domain.  The 
vitreous  domain  as  a  function  of  the  concentration  for 
B203-xLi20-yLi2S04  is  shown  in  Figure  10. 

The  results  on  ionic  conductivity  reported  here  are  obtained 
by  complex  impedance  measurements  in  the  frequency  range  5  Hz  to 
13  MHz  with  a  Hewlett-Packard  4192  ALF  impedance  analyzer.  The 
samples  used  are  rectangular  slabs  of  dimensions  4x2. 5x1  mm. 
Platinum  electrodes  of  circular  configuration  are  vapor  deposited 
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Figure  8.  Variation  of  the  Li+  conductivity  as  a  function  of 
the  reciprocal  temperature  for  the  borate  glasses 
B2°3”xL*2°"yL*2S04  ^or  different  concentrations  x  of 
L^O  and  y  of  L^SO^  according  to  reference  10.  The 
same  variation  for  B2O3-O.7  Li2O-0.7  LiCl  is 
represented  by  the  broken  line  (A.  Levasseur,  J.  C. 
Brethous,  J.  M.  Rfeau  and  P.  Hagenmuller,  Mat.  Res. 


riation  of  the  Li+  conductivity  at 
Oj-xL^O-yL^SO^  as  a  function  of 


on  the  two  sides  of  the  slab  as  shown  in  Figure  11.  The  applied 
voltage  along  the  z  direction  is  500  mV. 

1.2.2  Transport  properties  of  B2°3”0  .SL^O-yL^SO^ 

The  logarithm  of  the  conductivity  as  a  function  of  the 
reciprocal  temperature  for  the  different  samples  investigated  is 
plotted  in  Figure  12.  At  the  lower  temperature  side  the 
conductivity  follows  the  law 

-  -  £  «p  c-  J!f)  • 

The  deviations  from  this  law  observed  at  high  temperature  may  be 
due  to  the  difficulty  of  having  a  correct  reading  of  the 
temperature  at  high  temperatures. 

The  best  conductivity  measured  is  that  of  the  samples  with  a 

-2  -1  -1 

Li2S04  concentration  of  y  =  0.15:  <Jroax  ~  10  Q  cm  at 
30 0°C .  The  activation  energy  for  this  sample  is  Eg  =  0.6  eV. 
1.2.3.  Theoretical  Analysis 
1.2. 3.1.  Introduction 

In  the  system  that  we  are  dealing  with,  a  salt  such  as 
L^SOj  is  dissolved  in  a  solid  matrix  consisting  of  a  lithium 
borate  glass.  The  salt  dissociates  giving  Li+  ions  and 
LiS04  ions.  The  Li+  ions  are  highly  mobile  and  significantly 
increase  the  electrical  conductivity  of  the  system.  The 
situation  is  quite  analogous  to  an  electrolytic  solution  where  a 
salt  is  dissolved  in  a  liquid  solvent  such  as  water.  There  is 
one  significant  difference,  however,  and  that  is  that  water  has  a 
dielectric  constant  of  ~  80  whereas  borate  glasses  have  a 
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Figure  11.  Sample  configuration  for  the  complex  impedance 
measurements  on  B203“xt*i20-yLi2 SO^  showing  the 
position  of  the  platinum  electrodes. 
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LOG  o-T 


1.40  1.62  1.84  2.06  2.26 


1000/T  ( K “ 1 ) 


Figure  12.  Variation  of  log  a T  as  a  function  of  the  reciprocal 
absolute  temperature  for  different  concentrations  y 
of  L^SOj  in  82O3— 0.5  Li20— yLi2 SO^  • 


2 


[•] 


dielectric  constant  of  ~  o.  This  means  that  at  room  temperature 
a  typical  salt  such  as  LijSO^  will  be  nearly  fully  dissociated  in 
aqueous  solution  but  will  be  only  weakly  dissociated  in  the 
borate  glass.  Furthermore,  in  the  latter  case,  the  concentration 
of  lithium  ions  will  increase  rapidly  as  the  temperature  is 
increased  above  room  temperature.  Consequently,  the  electrical 
conductivity  of  the  doped  borate  glasses  should  increase  rapidly 
with  increasing  temperature,  as  seen  in  the  experimental  data 
discussed  in  Section  1.2  of  this  report.  In  the  following  sec¬ 
tions,  we  shall  develop  a  theory  for  the  temperature  dependence 
of  the  free  lithium  ion  concentration  and  derive  an  expression 
for  the  activation  energy  characterizing  the  conductivity. 

1.2. 3.2.  Calculation  of  mobile  Li*  ion  concentration 

Consider  an  ensemble  of  molecules  dispersed  in  the 

borate  glass  that  is  assumed  to  be  a  continuous  medium.  If  we 
regard  the  molecules  as  isolated  and  assume  that  ‘only  one  Li+  ion 
can  be  ionized  per  molecule  (these  assumptions  being  valid  in  the 
low  concentration  limit  at  temperatures  which  are  not  too  high) 
the  energy  of  an  Li+  ion  is  given  by 

2 

E(P,r>  •  +  v(t)  (1) 

where  V(r)  is  the  Coulomb  potential  energy  of  the  Li+  ion  in  the 
field  of  the  LiSO^  given  by 

2 
e 

c  r 


V(r) 


(2) 


p  is  the  momentum  of  the  Li+  ion,  r  is  its  distance  from  the 
negative  ion,  M  is  the  ion  mass  and  e  is  the  static  dielectric 
constant  of  the  glass.  The  number  of  ions  with  energy  E  is  given 
by  the  Boltzmann  distribution 

n  (E)  =  A  exp  {-  E/P;r)  }  (3) 

KB  1 

where  T  is  the  absolute  temperature  and  A  is  a  normalization 
constant  determined  by  the  total  number  N  of  ions  according  to 
the  expression 

M  =  A  /  d3p  J  d3r  exp  {-  }  (4) 

We  wish  to  calculate  the  number  NB  of  bound  Li+  ions  and  the 
number  of  free  Li+  ions  which  are  given  by 

N0  =  A  /  d3p  /  d3r  exp  {-  }  (5) 

for  (E  , p, r ) <0  and 

Np  =  A  /  d3 p  /  d3 r  exp  {-  }  (6) 

for  E(p,r)>0. 

From  Eqs .  (1)  and  (2)  we  can  see  that  the  conditions 
E(p,r)<0  and  E(p,r)>0  correspond  to  p2/2M<e2/er  and  p2/2M>e2/er 
respectively.  Introducing  spherical  coordinates  in  both  the  p 
and  r  spaces,  we  can  carry  out  the  integrals  over  the  angle 
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variables  and  obtain  the  expression 


Ng  =  16n1//2A(2nMknT)  3/2\3[-  ^.{  (  5/2 


-  i°)5/2\  ♦  I 


1/2  X/r0 

i —  J  y 

\/R 


-4 


exp(y)Erf (y1/2)dy 


(7) 


2 

where  X  =  e  /ek^T,  Erf(x)  is  the  error  function,  rQ  is  the 
distance  of  closest  approach  of  an  Li+  ion  to  an  LiSO^  ion,  and  R 
is  the  upper  limit  on  the  r  integration,  the  value  of  which  will 
be  discussed  shortly. 

We  can  evaluate  the  right-hand  side  of  Eq.  (6)  in  a  similar 
fashion  to  obtain  an  expression  for  Np.  This  gives 


=  16u1/2A(2nMkBT)  3/2X3[^Ur)  5/2 


(  r0  -j  5/2 1  n1/2  X{r0  -4  ,  . 

irH  /  +  -a —  1  y  e*p(y) 

q  \/R 


x  Erfc(yL/2)dy] 


(8) 


where  Erfc(x)  is  the  complement  of  the  error  function  given  by 
Erfc(x)  =  1  -  Erf(x).  It  can  be  shown  quite  easily  from  Eqs. 
(4),  (7)  and  (8)  that 


Nb  +  Np  =  N 


(9) 


23 


We  now  turn  to  the  discussion  of  the  upper  limit  R  and  the 

lower  limit  rQ  of  the  r  integration.  Suppose  that  initially  an 

Li+  ion  is  close  to  a  given  LiSO^  ion.  As  the  Li+  ion  moves 

further  away  from  the  LiSO^  ion  it  will  eventually  come  closer  to 

a  second  LiSO~  ion  and  will  then  move  predominantly  under  the 

influence  of  this  second  ion.  The  upper  limit  R  should  therefore 

be  approximately  given  by  R  =  r^/2  where  r^  is  the  mean  nearest- 

neighbor  distance  between  LiSO^  ions.  In  our  picture  the 

concentrations  of  the  free  Li+  ions  and  the  LiSOT  ions  have  the 

4 

same  value  nF  in  terms  of  which  we  can  express  r^  using  the 
relation 


=  nF~ 1/3  (10) 

Our  subsequent  analysis  of  the  experimental  data  will  show 
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that  n„  “  10  cm  and  so  r.  *  5.0x10  cm. 

The  lower  limit  rQ  can  be  estimated  from  Pauling's 
values*11^  for  the  S-0  distance  and  the  radii  for  0  and  Li+. 
Adding  these  values  gives  the  result  Tq  =  2.70A. 

We  also  need  to  make  an  estimate  of  the  quantity  \  which 
appears  in  Eqs.  (7)  and  (8) .  For  the  experimental  temperatures 
of  interest  (T  *  600K)  and  for  e  =  6.5,  a  typical  value  for 
lithium  borate  glasses,  we  find  that  \  "  4.3  x  10”^  cm.  It 
should  be  noted  that  \  and  r^  are  of  the  same  order  of  magnitude 
and  that  they  are  an  order  of  magnitude  larger  than  rg. 

Using  these  appraisals,  i.e.,  that  \  *  r^  >>  rQ,  we  can  make 
approximate  evaluations  of  the  integrals  in  Eqs.  (7)  and  (8)  to 
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obtain  the  results 


Nb  -  4nA(2nMkBT)  3/2X3(^)  4exp(i-)  (11) 

and 

Np  =  4itA(2nMkBT)3/2\3(ri/2\)3/it1/2  .  (12) 

Eliminating  A  from  Eqs .  (11)  and  (12)  and  using  Eq.  (10)  give  the 
relation 


N 


F 


-\/r 

NBn\e 


(13) 


where 


n\  ■ 


o 


(14) 


Dividing  both  sides  of  Eq.  (13)  by  the  volume  V  of  the  sample,  we 
get 


K(T) 


(15) 


where  Np  =  Vp/V  and  nB  =  NQ/V.  Equation  (15)  is  simply  the  mass 
action  law  for  the  reaction 


Li2S04  ♦  Li +  +  LiS0~ 


25 


with  K (T )  the  equilibrium  constant. 

At  this  point  we  must  take  into  account  the  fact  that  the 
conductivity  is  non-zero  when  y  =  0  due  to  a  small  concentration 
of  free  Li+  ions  arising  from  the  dissociation  of  Li20.  When  y 
becomes  different  from  zero,  there  will  still  be  some  concentra¬ 
tion  n£,  of  free  Li+  ions  from  the  Li20  which  must  be  included  in 
the  total  free  Li+  ion  concentration  nFT.  We  proceed  by  writing 
down  expressions  for  the  numbers  of  bound  Li+  ions  and  free 
Li+  ions  N£,  that  are  analogous  to  Eqs.  (7)  and  (8).  The  results 
are 

N'b  =  16h1/2A*  (2*MkBT) 3/2\3 {-(1/5) [ (£) 5/2-(^° ) 5/2 ] 

Vr; 

+  (*1/2/4)  /  y"*4eyErf  ( y 1/2 )  d y }  (16) 

\/R 


and 


N'  =  16it1/2A'  (27iMkBT)  3/2X3{  (1/5)  [  (^-)5/2-(^)5/2] 

F 

X/r^ 

+  (itL/2/4)  /  y“4eyErfc(y1/2)dy}  (17) 

X/  R 

where  r^  is  the  distance  of  closest  approach  of  a  Li+  ion  to  a 
LiO“  ion.  A'  is  determined  from  the  relation  +  N^,  =  N' ,  and  N' 
is  the  total  concentration  of  Li+  ions  that  can  be  ionized  from 
Li20  according  to  the  reaction  L^O  >  Li+  +  LiO-.  From  Pauling's 
radii  we  obtain  the  value  of  2.0A  for  r^. 

We  can  obtain  an  approximate  expression  for  N£,  by  using  the 
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same  procedure  that  led  to  Eq.  (12) .  The  quantity  r^  is  now  to 
be  related  to  the  mean  nearest  neighbor  separation  between  any 
pair  of  negative  ions,  either  LiO“  or  LiSO”,  and  is  taken  to  be 

r  i  =(np  +  np-1/3  .  (18) 

Taking  R  =  r ^/2  and  using  the  appraisals  \  -  r^  >>  r^,  we  find 
that 

Mp  3  4nA'  (2nMk0T)  3/2  (ri/2)  3/u1/2  .  (19) 

The  values  of  r^  given  by  Eq.  (18)  should  also  be  used  in  Eq. 

(12)  . 

We  now  proceed  in  a  fashion  analogous  to  that  which  led  to 
Eq.  (15).  The  mass  action  laws  for  the  two  reactions 

Li2S04  *■  Li+  +  LiSO^ 

Li20  *  Li  +  +  LiO* 

take  the  form 


'B 


i 


nFTno 


-  K(T) 


‘B 


"FT% 


1  ^ro 

nT  e 


K'  (T) 


(20) 


(21) 


where  n?T  *  np  +  np,  nQ  is  the  concentration  of  LiSO^  ions,  n^  is 
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the  concentration  of  LiO“  ions,  n£  is  the  concentration  of  Li+ 
ions  bound  to  LiO”  ions,  and 


n 


X 


X 

1 7? 


8n 


<r;> 


(22) 


Let  us  introduce  the  quantities  n  =  nQ  +  nB  and  n'  =  n^  +  n^ 
which  represent,  respectively,  the  total  concentrations  of  Li2S04 
and  Li20.  If  we  eliminate  nB  and  n^  from  the  expressions  for  n 
and  n'  using  Eqs .  (21)  and  (22)  and  rearrange  the  results,  we  get 


"o  i+nF^K(T) 
no  =  l+np^K' (TV 

From  conservation  of  lithium,  we  can  write 


nFT 


n. 


+  n 


(25) 


and  eliminating  nQ  and  n^  with  the  aid  of  Eqs.  (24)  and  (25),  we 
obtain 


nFT  “  r+'n'F"K rry  +  r+n^rnrr  *  (2 

Since  we  expect  that  the  quantities  nF7/n  and  nFT/n'  are  small, 
we  can  assume  that  nFTK(T)>>l  and  nFTK'(T)>>l.  Equation  (26) 
then  simplifies  to 
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FT  K  (T )  K 1  (T)  *  U/ 

Substituting  the  expressions  for  K(T)  and  K'(T)  from  Sqs.  (20) 
and  (21)  into  Eq .  (27),  and  taking  the  square  root  of  both  sides 

leads  to  the  result 


"ft  =  tyn\nTe 


-2Ead/kBT 


+  xn\nTe 


-2E»d/kBT]l/2 


(28) 


where  x  =  n'/nT, 


and 


p*  -  e 
Ead  _  7TT 
o 


(29) 


(30) 


The  quantities  Ead  and  are  the  limiting  activation  energies 

for  dissociation  associated  with  the  cases  x  =  0  and  y  =  0, 
respectively. 

Using  the  expression  a  =  nFTevi,  the  conductivity  can  be 
written  as 


a  =  eian 


1/2 


[yn 


-2e»<i/|'bt  + 

e  + 


xn 


»  e”2Ead/kBT j 1/2 


(31) 


Using  the  Einstein  relation  connecting  the  mobility  n  and  the 
diffusion  coefficient  D  and  recognizing  that  both  quantities  may 
have  an  activation  energy  Egh  associated  with  activated  hopping, 
we  can  write 
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e 


(32) 


eD 

-  T^T 


eDo 


“Eaf/kBT 


Substituting  Eq.  (32)  into  Eq .  (31)  and  using  definitions  of  n^ 
and  n^,  we  see  that  a  plot  of  ln(aT  '  )  versus  1/T  for  x  =  0  or 
y  =  0  should  yield  a  straight  line.  Even  when  neither  x  nor  y  is 
zero,  this  plot  should  still  be  very  close  to  a  straight  line. 

The  slope  of  the  line  multiplied  by  -kB  gives  the  effective 
activation  energy  for  the  conductivity,  Ea. 

The  system  that  we  are  considering  can  be  compared  to  that 
of  an  alkali  halide  such  as  NaCl  doped  with  CaCl2  or  CdCl2,  which 
has  been  discussed  previously . ^ ^  The  Ca++  or  Cd++  ions 
substitute  for  Na+  ions  in  the  lattice,  and  to  preserve 
electroneutrality  positive  ion  vacancies  must  be  introduced 
simultaneously  with  a  concentration  equal  to  that  of  Ca++  or 
Cd++.  Positive  ion  vacancies  behave  as  negative  charges  which, 
when  dissociated  from  the  Ca++  or  Cd++  ions,  are  mobile  and  can 
contribute  to  the  electrical  conductivity.  The  overall 
activation  energy  for  the  conductivity  is  the  sum  of  a  part 
connected  with  the  dissociation  of  a  positive  ion  vacancy  from  a 
divalent  cation  and  a  part  connected  with  the  diffusion  of  a 
positive  ion  vacancy  under  the  influence  of  an  applied  electric 
field.  In  our  system  the  mobile  lithium  ions  are  the  analogues 
of  the  positive  vacancies.  In  the  standard  treatment  a  positive 
ion  vacancy  is  regarded  as  completely  free  unless  it  occupies  one 
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of  the  twelve  nearest-neighbor  neighbor  sites  of  a  divalent 
cation.  We  make  no  such  assumption  in  our  treatment  and  take 
into  account  the  Coulomb  interaction  at  large  distances  of  a  Li+ 

ion  from  a  LiO"  or  LiSOT  ion. 

4 

1.2. 3. 3.  Numerical  Results  and  Comparison  with  Experiment 

Numerical  calculations  of  the  conductivity  as  a  function  of 
temperature  and  composition  were  carried  out  in  two  ways:  (i) 
the  conductivity  was  calculated  from  the  relation  a  =  (n^+npeii 
using  the  exact  expressions  contained  in  Eqs .  (7)  and  (8)  and  in 
Eqs .  (16)  and  (17)  to  determine  nF  and  n£,  respectively;  and  (ii) 

the  conductivity  was  calculated  from  Eq.  (31)  which  is  based  on 
the  approximations  \  =  ri  >>  rQ,  r^.  In  both  cases  the  parameter 

no  o 

nT  was  taken  to  be  1.0x10^  cm-J  and  the  parameters  DQ,  Egh,  and 
e  were  varied  to  give  the  best  least  squares  fit  to  the 
experimental  data.  The  mobility  activation  energy  Eatl  was 
determined  from  the  y=0  data  (only  Li20  present)  to  be  0.38  eV 
and  was  kept  at  this  same  value  for  the  ytO  data  (Li2S04 
present) .  The  parameters  DQ  and  e  were  determined  for  each  value 
of  y . 

The  results  are  the  following. 

Case  (i) 

3/2 

The  quantity  ln(oT  )  is  plotted  against  the  reciprocal  of 
the  absolute  temperature  in  Figure  13  for  each  of  the  composi¬ 
tions  y  =  0.0,  0.05,  0.10,  and  0.15.  The  theoretical  curves  are 
straight  lines  which  fit  the  experimental  data  quite  well  except 
for  some  deviations  at  the  two  lowest  temperatures.  From  the 
values  of  e  determined  by  the  fits  to  the  data,  the  values  of  the 
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Figure  13.  Variation  of  log  (oT  '  )  as  a  function  of  the 
reciprocal  absolute  temperature  for  different 
concentrations  y  of  Li2S04  in  B2O3-O . SL^O-yL^SO^ . 
The  theoretical  results  are  represented  by  the  solid 
curves . 
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CONDUCTIVITY  cr  (10'°^'' cm 


CONCENTRATION  OF  Li2S04  y 


Figure  14a.  Variation  of  the  conductivity  a  as  a  function  of 
Li2S04  concentration  y  for  T  »  469  K  and  489  K. 
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CONCENTRATION  OF  Li2S04  y 

Figure  14b.  Variation  of  the  conductivity  a  as  a  function  of 
Li2S04  concentration  y  for  T  ■  539  K  and  569  K. 
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CONCENTRATION  OF  Li2S04  y 

Figure  14c.  Variation  of  the  conductivity  a  as  a  function  of 
L^SOj  concentration  y  for  T  ■  619  K  and  659  K. 
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dissociation  activation  energies  Ead  and  were  calculated 

using  Eqs.  (29)  and  (30) .  The  overall  activation  energy  Eg  was 
determined  from  the  slopes  of  the  theoretical  curves  in  Figure 
13.  The  results  for  DQ,  e,  Ea(g,  Ead '  an<^  theoretical  and 

experimental  values  of  EaA  are  presented  in  Table  2.  We  see  that 
the  four  quantities  DQ,  e,  and  Egd,  and  E^d  show  only  a  weak 
dependence  on  the  composition  variable  y,  and  their  magnitudes 
are  physically  reasonable.  That  of  DQ  conforms  roughly  to  the 
relation ^12^  Dq  =  v+r^/6  where  v+  and  r^  are  the  hopping 
frequency  and  hopping  distance,  respectively.  Taking  0.2  cm2 
sec-1  as  a  representative  value  of  DQ  from  Table  2  and  taking^2) 
v+  =  10  sec  ,  we  find  that  r^  *  11  A,  a  not  unreasonable 
value. 


The  conductivity  a  is  plotted  against  y  in  Figure  14  for  the 
temperatures  469K,  489K,  539K,  569K,  619K,  and  659K.  As  in 
Figure  13  there  are  appreciable  deviations  between  theory  and 
experiment  for  the  two  lowest  temperatures,  but  for  the  four 
highest  temperatures  the  agreement  is  good.  One  can,  of  course, 
improve  the  agreement  for  the  two  lowest  temperatures  by  making  e 
temperature  dependent. 

Case  (ii) 

When  the  conductivity  a  is  calculated  using  Eq.  (31) ,  the 
3 /2 

plots  of  ln(aT  '  )  versus  1/T  fall  within  one  percent  of  those 
given  in  Figure  13  for  case  (i)  for  each  value  of  y.  Plots  of  a 
versus  y  also  fall  within  one  percent  of  those  given  in  Figure  14 
for  case  (i).  These  results  indicate  that  the  approximations 
underlying  Eq.  (31)  are  well  satisfied. 
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TABLE  2 


CONDUCTIVITY  PARAMETERS  FOR  B203-xL i 20-yL i 2S04 


y 

Do 

e 

Ead 

Ead ' 

Ea 

(theor . ) 

Ea 

(exp.) 

0.00 

0.15 

6.6 

0.55 

0.93 

0.93 

0.05 

0.26 

6.0 

0.44 

0.60 

0.83 

0.90 

0.10 

0.21 

6.5 

0.41 

0.55 

0.79 

0.83 

0.15 

0.21 

6.6 

0.40 

0.55 

0.78 

0.83 
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The  theory which  we  have  developed  for  the  dissociation 
of  L^SO^  in  lithium  borate  glasses  *3203-1^20  gives  a  good 
representation  of  the  dependence  of  the  conductivity  on  the 
Li 2 SO 4  content  and  on  the  temperature.  By  comparing  the 
experimentally  determined  conductivity  with  the  theoretical 
results  one  obtains  information  concerning  the  activation  energy 
associated  with  the  mobility,  the  diffusion  constant  and  the 
dielectric  constant. 

1.2. 3. 4.  Coulomb  Interaction  between  Free  Li+  Ions 

The  effect  of  the  Coulomb  interaction  between  free  Li+  ions 
can  be  taken  into  account  by  using  the  Debye-Huckel  theory  in  a 
form  applied  by  Lidiard^^  to  a  similar  problem.  What  has  to  be 
done  is  to  replace  the  dissociation  constants  K (T)  and  K' (T)  by 
Kdf(T)  and  K^fT)  given  by 

KDH(T'  ■  HT  -  ^<2M3/2(4.nFT)1/2  (33) 

K  O  O 

KlVT>  ■  KT  explrr  -  £[<2M3/2(4«nFT)1/2  Tfijr])  (34) 

A.  O  O 

2 

where  <  -  8unFTX..  Since  the  new  equilibrium  constants  involve 

nFT,  the  mass  action  law  equations  are  now  transcendental 
equations  for  the  determination  of  nFT  and  must  be  solved  self- 
consistently . 

If  we  introduce  the  quantities  n  *  nQ  +  nB  and  n'  =  n^  +  n^ 
which  are,  respectively,  the  total  concentrations  of  Li2S04  and 
Li20,  we  can  rewrite  Eqs .  (23)  and  (24)  in  the  form 
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°  1  +  nFTKDH  ^ 


n 


(35) 


n* 


no  - 


(36) 


Using  conservation  of  lithium  ions,  we  have 


'FT 


=  nrt  +  n' 
o  o 


1+nFTKDH(T^ 


+ 


1+nFTKDH (T^ 


(37) 


Equation  (37)  is  the  equation  that  must  be  solved  self 
consistently  to  give  the  total  free  Li+  ion  concentration  np^,. 

Once  nFT  is  found,  we  can  determine  the  conductivity  a  using 
the  relation 


®  ~  en  n  pip 


(38) 


where  p.  is  the  mobility.  In  the  presence  of  the  Coulomb 
interaction  between  the  free  Li+  ions,  the  mobility  pn  for  the 
noninteracting  case  must  be  multiplied  by  a  factor  C  given  by 


C 


nf  +  n'  f ' 
n+n'  n+n' 


(39) 


where 


f  =  1 


_ _ 

3(/2+l)(l+Kr0)(/2+Kr0) 


(40a) 
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f '  =  1 - = - — - — -  .  (40b) 

3(/2  +  l)(l+Kr^)(/2+K:r^) 

The  expression  for  the  conductivity  then  takes  the  form 

a  =  ennCnFT  .  (41) 

We  have  used  Eqs .  (37)  and  (41)  to  calculate  the 

conductivity  for  the  system  B2O3-O . 5Li 20-yLi 2SO4  as  a  function  of 

y  and  the  temperature  T.  The  values  of  the  parameters  e,  rQ,  and 

r'  used  here  are:  e  =  6.5,  r  =  2.70A,  r'  =  2.0A.  The  results 
o  '00 

are  presented  in  Figures  15,  16,  and  17  for  the  various 
compositions  and  temperatures  considered.  At  the  two  lowest 
temperatures,  there  is  a  small,  but  distinct  lowering  of  the 
conductivity  that  arises  from  the  Coulomb  interaction  between  the 
free  Li+  ions.  A  given  ion  is  slowed  up  by  its  interaction  with 
its  ion  atmosphere.  As  the  temperature  increases,  the  lowering 
of  the  conductivity  diminishes  and  is  hardly  evident  at  the  two 
highest  temperatures.  This  behavior  can  be  understood  on  the 
basis  that  an  increase  in  temperature  increases  the  random 
thermal  motion  of  the  free  Li+  ions  and  reduces  the  correlations 
that  arise  from  the  Coulomb  interaction  between  them.  These 
results  are  consistent  with  the  picture  that  the  lithium  borate 
glasses  behave  as  weak  electrolytes  with  only  a  small  fraction  of 
the  L^O  and  L^SO^  molecules  ionized  at  the  temperatures 
employed  in  the  experiments. 
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CONDUCTIVITY  a  (I03&  'cm 


0.00  0.04  0.08  0.12  0.16 

CONCENTRATION  OF  Li2S04  y 

Figure  15.  Variation  of  the  conductivity  a  as  a  function  of 

Li2S04  concentration  y  for  T  ■  469  K  and  489  K.  The 
solid  and  dashed  lines  are  the  cases  with  and  without 
inclusion  of  the  interaction  between  free  Li+  ions, 
respectively. 
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CONDUCTIVITY  a  (IO’*a"'cm 
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CONCENTRATION  OF  Li2S04  y 

Figure  16.  Variation  of  the  conductivity  a  as  a  function  of 

Li2S04  concentration  y  for  T  »  539  K  and  569  K.  The 
solid  and  dashed  lines  are  the  cases  with  and  without 
inclusion  of  the  interaction  between  free  Li+  ions, 
respectively. 
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CONCENTRATION  OF  Li2S04  y 


Figure  17.  Variation  of  the  conductivity  a  as  a  function  of 

Li2S04  concentration  y  for  T  »  619  K  and  659  K.  The 
cases  with  and  without  inclusion  of  the  interaction 
between  free  Li+  ions  are  indistinguishable. 
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1.3. 


Lithium  Ion  Vibration  in  Borate  Glasses 


In  the  discussion  of  the  dynamics  of  the  Li+  ions  we  shall 
refer  to  the  borate  glass  structure  shown  m  Figure  1. 

Reflectivity  measurements  were  performed  with  a  Bruker 
Infrared  Fourier  Transform  Spectrometer  (IFS)  at  room 
temperature.  The  spectral  range  analyzed  is  50  to  2000  cm"^,  the 
sample  used  on  parallelepipeds  of  dimensions  8x8x4  mm. 

Infrared  reflectivity  spectra  for  92°3“0  •  TL^O-yL^SC^  at 
T=300K  are  shown  in  Figure  18.  In  the  spectrum  corresponding  to 
the  composition  x=0.7,  y=0  three  main  features  are  identified  at 
1400  cm--*-,  1000  cm"^,  and  700  cm“^  associated  with  the  infrared 
active  vibrational  modes  of  803,  BO4  and  B-O-B,  respectively. 
Addition  of  the  dopant  L^SC^  introduces  two  new  bands  at  600 
cm-1  and  1100  cm""^  associated  with  the  and  modes  of  the 
sulfate  ion  SO~“ : 
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40 


Figure  18.  Infrared  reflectivity  spectra  for 
B2O3-O ,7Li20-yLi2S04  at  T  *  300K. 
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The  band  which  corresponds  to  the  Li  vibration  at  400  cm“^ 
is  covered  by  a  broad  band  probably  due  to  the  molecular 
reorientation  in  the  system  (Debye  relaxation) . 

In  order  to  analyze  the  reflectivity  spectra  it  is  necessary 
to  perform  a  Kramer s-Kronig  transformation  to  get  the  optical 
constants  associated  with  each  of  the  vibrational  modes.  In 
Figure  19  is  shown  the  imaginary  part  of  the  dielectric  constant 
associated  with  the  transverse  optical  phonon  (TO)  mode.  The 
vibrational  frequency  of  a  Li+  ion  is  situated  at  380  cm“l.  This 
vibration  concerns  essentially  the  Li+  -0  bonds  of  the  vitreous 
matr ix . 

In  Figure  20  is  presented  the  frequency  dependent 
conductivity  which  shows  a  resonance  at  a  frequency  close  to  380 
cm-1.  The  shapes  of  the  spectra  for  doped  and  undoped  samples 
are  similar  which  is  consistent  with  the  fact  that  B2C>3-Li20  is 
also  a  Li+  conductor.  This  result  also  shows  that  there  is 
practically  no  difference  in  the  structure  of  B2O3-0.7Li2O  3 

b2o3-0 «7Li2O-0 .4Li2S04  which  indicates  that  for  concentrations 
y<0.4  the  dominant  Li  sites  are  the  same. 

2 •  LITHIUM  intercalation  in  mse 

2.1.  Characteristics  of  the  intercalation  Process 

9 

2.2.1.  Introduction 

Lithium  intercalation  in  InSe  single  crystals  has  been 
performed for  the  first  time  in  the  Laboratoire  de  Physique 
des  Solides  of  the  Universite  Pierre  et  Marie  Curie,  Paris. 
Chemical  insertion  of  Li  in  InSe  is  obtained  by  insertion  of  a 
single  crystal  of  InSe  in  a  1.6  M  solution  of  n-butyll ithi urn  in 
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DIELECTRIC  CONSTANT 


Figure  19.  Imaginary  part  of  the  dielectric  constant  for 
B203-0.7Li20-yLi2S04  at  T  *  300  K. 
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Figure  20.  Frequency  dependent  conductivil 

L20-yLi2S04  . 


B203-0 .7L i 


hexane.  Spontaneous  intercalation  is  produced  under  these 
conditions  according  to  the  reaction. 

n  -  C4HgLi  +  InSe  ■+  LixInSe  +  j  n-C8H18* 

The  reaction  lasts  for  a  long  time  under  an  inert  atmosphere. 

The  following  steps  should  be  distinguished  in  this 
react i on . 

1)  pre-equilibrium  in  the  solution  generating  the  active 
soecies:  monomer,  dimer  and  trimer  of  n-C^HgLi 

ii)  surface  adsorption  and  activation  of  the  adsorbant 
(n-C4H9Li) 

iii)  electron  transfer  e  -*■  Li  bound  to  the  solid 

iv)  diffusion  of  Li  into  the  lattice 

v)  ionization  of  Li 

vi)  d imer isat ion ,  diffusion  of  the  alkane  in  the  solution 
and  production  of  octane. 

The  positions  of  the  Li  ions  in  the  Van  der  Waals  gaps  of 
InSe  are  shown  in  Figure  21: 

No  indication  of  staging  has  been  observed  in  this  process 
of  intercalation,  but  a  saturation  as  a  function  of  time  has  been 
observed.  When  lithium  is  intercalated  into  the  layer  compound 
InSe,  electrical  conductivity  measurements  indicate  that  the  rate 
of  intercalation  is  large  at  first,  then  decreases,  and  finally 
approaches  zero.  This  behavior  indicates  that  there  is  a  maximum 
or  saturation  value  of  the  lithium  concentration  that  can  be 
achieved.  A  physical  picture  of  this  saturation  behavior  is  as 
follows.  The  lithium  atoms  enter  the  Van  der  Waals  gap  between 
adjacent  InSe  layers.  The  2S  electrons  of  the  lithium  atoms 
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Positions  of  Li+  io 
Waals  gap. 


ionize  and  transfer  to  a  two-dimensional  potential  well  in  one  of 
the  adjacent  InSe  layers.  As  more  lithium  atoms  are 
intercalated,  the  Fermi  energy  Ep  of  the  electrons  in  the 
ootential  well  rises  as  a  consequence  of  the  Fermi-Dirac 
statistics  obeyed  by  the  electrons.  When  the  Fermi  energy 
becomes  equal  to  the  energy  of  a  2S  electron  in  a  lithium  atom, 
the  transfer  of  electrons  from  the  atoms  to  the  potential  well 
will  cease  and  the  electrical  conductivity  due  to  the  electrons 
in  the  potential  well  will  saturate. 

2.1.2.  Time  Evolution  of  Intercalation 

The  intercalation  of  Li  into  InSe  is  achieved  by  immersing  a 
single  crystal  of  InSe  into  a  1.6  M  solution  off  n-butyll i thium 
in  hexane.  Experimental  measurements  of.  the  resistivity  of  the 
InSe  during  intercalation  have  been  made  by  C.  Julien  and 
coworkers  in  the  laboratory  of  Professor  M.  Balkanski  in  Paris. 
The  results  for  the  resistivity  as  a  function  of  time  of 
intercalation  are  shown  in  Figure  22.  It  can  be  seen  that  the 
resistivity  initially  drops  quite  rapidly  and  then  levels  off. 

One  might  surmise  that  saturation  is  being  approached  at  the 
largest  times.  However,  this  does  not  seem  to  be  the  case.  One 
must  beware  of  the  logarithmic  scale  used  on  the  vertical  axis  in 
Figure  22.  If  one  plots  the  reciprocal  of  the  resistivity,  or 
conductivity,  versus  time  of  intercalation,  one  obtains  the  plot 
shown  in  Figure  23.  We  see  that  the  conductivity  is  very  nearly 
a  linear  function  of  time  over  the  interval  considered  and  that 
no  evidence  of  saturation  is  shown.  Viewed  in  this  light,  the 
experimental  data  are  consistent  with  the  picture  that  the 
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conduction  electron  concentration  in  the  InSe  is  increasing 
linearly  with  time  as  a  consequence  of  a  corresponding  linear 
increase  in  Li  concentration  with  time. 

2.1.3.  Analysis  of  Saturation  of  Electrical  Conductivity  During 
Intercalation 

Although  saturation  of  the  conductivity  has  not  yet  been 
achieved  experimentally,  it  should  be  attainable  after 
sufficiently  long  periods  of  intercalation.  In  our  physical 
picture,  the  intercalated  Li  atoms  enter  the  Van  der  Waals  gap 
between  adjacent  InSe  layers  and  lose  their  2S  electrons  to  one 
or  the  other  of  the  two-dimensional  potential  wells  in  the 
adjacent  InSe  layers.  As  more  Li  atoms  are  intercalated,  the 
Fermi  energy  Ep  of  the  electrons  in  the  potential  well  rises  as  a 
consequence  of  the  Fermi-Dirac  statistics  obeyed  by  the 
electrons.  When  the  Fermi  energy  becomes  equal  to  the  energy  of 
a  2S  electron  in  a  lithium  atom,  the  transfer  of  electrons  from 
the  atoms  to  the  potential  well  will  cease  and  the  electrical 
conductivity  due  to  the  electrons  in  the  potential  well  will 
saturate . 

We  have  carried  out  a  Thomas-Fermi  analysis  of  the  fraction 
of  electrons  transferred  as  a  function  of  the  well  depth.  Let 
the  width  of  the  Van  der  Waals  gap  be  aQ  and  the  periodicity 
length  perpendicular  to  the  planes  of  the  InSe  be  cQ.  The 
quantity  cQ-a0  is  therefore  the  width  of  the  potential  well. 
Within  the  well,  the  potential  has  the  constant  value  VQ.  The 
Li+  ions  are  assumed  to  be  uniformly  distributed  over  two  planes 
symmetrically  disposed  about  the  middle  of  the  Van  der  Waals  gap 
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and  a  distance  C+  apart. 

The  Thomas-Fermi  procedure  utilizes  an  energy  functional  of 
the  electron  number  density  which  consists  of  the  kinetic  energy 
and  electrostatic  potential  energy.  By  minimizing  the  energy 
function  with  respect  to  variations  in  number  density,  one 
obtains  a  relation  between  the  latter  and  the  electric 
potential.  Combining  this  relation  with  Poisson's  equation  and 
applying  the  appropriate  boundary  conditions  leads  to  an 
expression  for  the  number  density  as  a  function  of  position  in 
the  system.  The  fraction  of  the  electrons  located  within  the 
InSe  potential  well  can  then  be  readily  calculated,  and  the 
result  is  (details  of  the  calculation  can  be  found  in  R&D  status 
report  No.  3  for  contract  No.  F33615-85-K-2501) . 
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where  is  the  electron  charge  in  the  Van  der  Waals  gap,  Q2  is 

ao  h2 

that  in  the  InSe  layer,  «  =  — ,  aR  =  — =■  is  the  Bohr  radius,  e  is 

B  me 

★ 

the  mean  dielectric  constant  and  x  =  with  n  =  the  reduced 

co  « 

electron  mass  and  c  =  — . 

B 

The  fraction  f  increases  with  increasing  potential  well 
depth  VQ  until  a  critical  value  of  V0  is  reached  for  which  f  =  1 
and  all  of  the  electrons  have  been  transferred  from  the  lithium 
atoms  to  the  InSe  layer.  In  Figure  24  is  giv»n  a  plot  of  f  versus 
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Figure  24.  Charge  transfer  at  saturation  versus  potential  well 
depth  for  Li-InSe  intercalation. 
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V0/Pn*  We  see  that  f  increases  linearly  with  VQ .  The  critical 
value  of  VQ  is  found  to  be  Vcr^t  =  1.043  pQ.  The  minimum  value 
of  f  is  not  zero  because  the  Li+  ions  are  assumed  to  be  at  the 
interface  between  Van  der  Waals  gap  and  the  InSe  layer.  A  more 
realistic  calculation  is  currently  underway  in  which  the  Li+  ions 
are  moved  into  the  Van  der  Waals  gap  away  from  the  interface. 

2.2.  might  Binding  Calculation  of  the  Energy  Bands  in  InSe 

The  first  step  in  calculating  the  potential  energy  surface 
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The  activation  energy  is  the  energy  difference  between  a 
local  potential  energy  minimum  of  a  lithium  ion  and  the  saddle 
point  between  that  minimum  and  the  next  one.  It  is  therefore 
necessary  to  calculate  the  potential  energy  surface  for  the 
motion  of  a  lithium  ion  in  InSe  in  order  to  obtain  the  activation 
energy.  We  have  carried  out  such  a  calculation  using  the  Thomas- 
Fermi  approximation  in  which  the  potential  3>(r)  at  a  given  point 
in  the  crystal  is  related  to  the  electron  density  at  that  point 
by  the  relation 
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(43) 


®(r) 


where  Rj  is  the  position  vector  of  the  j-th  atomic  site  and  C  is 
the  normalization  constant.  The  atomic  orbitals  were  taken  to  be 
those  tabulated  by  Herman  and  Skillman^^,  and  the  tight  binding 
parameters  were  taken  to  be  those  of  Doni  et  al.(17). 

Using  the  Bloch  functions  for  ^  =  0,  the  electron  density 
p(r)  was  calculated  using  Eq.  (44).  The  potential  $(r)  was  then 
calculated  using  Eq.  (43).  The  potential  energy  V(r)  is  given  by 
V(r)  =  e$(r) .  A  contour  plot  of  the  potential  energy  as  a 
function  of  position  in  the  unit  cell  is  given  in  Figure  25  for 
the  plane  lying  in  the  center  of  the  Van  der  Waals  gap  between 
successive  InSe  layers.  The  positions  of  minimum  potential 
energy  are  indicated  by  squares,  and  the  saddle  point 
corresponding  to  the  minimum  activation  energy  for  ion  motion  is 
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Figure  25.  Potential  energy  contours  for  a  lithium  ion  moving  in 
the  plane  of  the  Van  der  Waals  gap  of  InSe.  The 
regions  of  high  potential  energy  are  selenium  atoms. 
□  :potential  minimum;  A;  saddle  point. 
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indicated  by  the  triangle.  The  contours  are  labeled  by  the 
energies  in  Rydbergs.  The  activation  energy  is  the  difference  in 
energy  of  the  sites  labeled  by  squares  and  triangles  and  is  found 
to  be  2.2  eV. 

We  feel  that  the  activation  energy  just  quoted  is  rather 
large  in  view  of  the  apparent  ease  with  which  intercalated 
lithium  atoms  move  through  the  inSe  lattice.  This  calculation 
was  repeated  by  taking  into  account  all  values  of  the  wave  vector 
t  in  the  first  Brillouin  zone.  A  new  contour  plot  is  obtained  as 
shown  in  Figure  26.  The  potential  energy  along  the  path  from  one 
minimum  to  the  next  through  the  saddle  point  is  shown  in  Figure 
27.  From  this  plot  the  activation  energy  can  be  seen  to  be  1.6 
eV,  which  is  significantly  smaller  than  the  value  of  2.2  eV 
obtained  using  only  the  &  =  0  Bloch  functions. 

The  1.6  eV  value  is  still  rather  large  in  view  of  the 
apparent  ease  with  which  intercalated  lithium  atoms  move  through 
the  InSe  lattice.  Since  we  have  up  to  now  neglected  the  2s 
electron  of  the  lithium  atom,  it  seems  appropriate  that  we 
include  it  in  our  calculations. 

3 .  FAST  ION  TRANSFER  BETWEEN  LITHIUM  BORATE  GLASS  AND  InSe  LAYER 
COMPOUND 

3.1.'  introduction 

In  the  solid  state  battery  concept  based  on  a  doped  lithium 
borate  glass  as  the  electrolyte  and  indium  selenide  intercalated 
with  lithium  as  the  cathode,  the  borate  glass  is  a  three- 
dimensional  fast  ion  conductor  whereas  the  indium  selenide  is  a 
two-dimensional  fast  ion  conductor  with  the  lithium  ions  moving 
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Figure  26.  The  new  potential  energy  contours  for  a  lithium  ion 


moving  in  the  plane  of  the  Van  der  Waals  gap  of  InSe 
when  all  values  of  the  wave  vector  $  in  the  first 
Brillouin  zone  are  taken  into  account.  The  regions 
of  high  potential  energy  are  selenium  atoms. 
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PATH  OISTANCE  (angstroms) 


Figure  27.  The  potential  energy  along  tfie  path  from  one  minimum 
to  the  next  through  the  saddle  point. 
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in  the  Van  der  Waals  gaps  between  the  InSe  layers.  The  ion 
current  must  therefore  make  a  transition  from  three-dimensional 
to  two-dimensional  behavior  as  the  Li+  ions  pass  from  the  glass 
into  the  InSe. 

In  Section  1.2,  we  presented  a  theory  of  the  electrical 
conductivity  of  a  lithium  borate  glass  doped  with  L^SO^.  The 
conductivity  as  a  function  of  temperature  and  composition  was 
calculated  and  found  to  be  in  good  agreement  with  experimental 
data.  In  this  chapter,  we  give  the  results  of  a  theoretical 
investigation  of  the  constant  current  contours  and  constant 
potential  contours  for  the  flow  of  positive  ions  from  the  glass 
into  the  InSe  under  the  influence  of  an  applied  voltage.  The 
effective  resistance  of  the  glass  as  a  function  of  the 
geometrical  parameters  is  also  calculated. 

3.2.  Theoretical  Development 

The  geometry  of  our  system  is  shown  in  Figure  28.  The  Van 
der  Waals  gaps  in  the  InSe  have  width  2c  and  the  distance  between 
successive  Van  der  Waals  gaps  is  2w.  A  slab  of  borate  glass  of 
width  2w  and  length  L  can  be  viewed  as  centered  against  each  Van 
der  Waals  gap.  A  potential  difference  A<|>  is  imposed  across  the 
borate  glass  as  shown.  We  require  that  the  Li+  ions  enter  the 
InSe  only  through  the  Van  der  Waals  gaps.  The  borate  glass  is 
assumed  to  have  a  conductivity  a. 

The  problem  we  address  is  to  calculate  the  current  density  J 
in  the  borate  glass  as  a  function  of  c,  w,  L  and  A<>.  We  proceed 
by  solving  Laplace's  equation 
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0 


(46) 


A  $ 


subject  to  appropriate  boundary  conditions.  If  we  assume  that 
the  system  extends  to  ±  “  in  the  y-direction,  we  can  exploit  the 
periodicity  of  the  system  by  utilizing  Fourier  series. 

The  boundary  conditions  which  the  potential  4>(x,y)  satisfies 
are  the  following: 


4>(0,y)  = 

1  ^a 

(47a) 

0(L,y)  . 

=  ^b 

(47b) 

9<t>  (x,y) 

,-v° 

x=L  ~  0 

for  | y |  <  c 

(47c) 

dx 

for  c  <  |  y  j  <  w 

2»<t>  (x,y) 

=  0  for 
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0  <  x  <  L 

(47d) 

where  JQ  is 
InSe  at  the 
solution  of 
in  the  form 


the  current  density  across  a  Van  der  Waals  gap  of 

interface  with  the  lithium  borate  glass.  The 

Eq.  (46)  which  satisfies  Eqs.  (47)  can  be  expressed 
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The  potential  difference  A$  *  $a  -  -  <|>(L,o)  is  given  by 
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Equation  (49)  establishes  the  relation  between  the  current 
density  JQ  and  the  applied  potential  difference  A 4> . 

The  x-  and  y-components  of  the  current  density  J  are  given 
by 


J  =  -  a 


S7 


J  [§  +  -  l  - 
oLW  n 


„  .  nitx  nit 

cosh  -  cos 


w 


w 


sin 


nnc 
w  ■ 


n-1 


cos 


nwL 

w 


J  =  -  a 


a  <t> 

3y 


=  J  [-  -  l  — 
o  L  n  ‘'in 
n=  1 


eiBV,  nitx  _j_  miy  nnc 

sinh  -  sm  — sin  - 

w  w  w  j 


cosh 


nnL 

W 


(50a) 


(50b) 


3.3.  Numerical  Results 

Numerical  calculations  of  Jx  and  Jy  as  functions  of  x  and  y 
have  been  carried  out  for  various  values  of  the  geometrical 
parameters,  c,  w,  and  L.  Plots  of  the  contours  of  constant  Jx 
and  Jy  (in  units  of  JQ)  are  presented  in  Figures  29  and  30  for 
the  case  c  =  0.2w  and  L  =  4w.  The  channeling  of  the  ion  current 
into  the  Van  der  Waals  gap  is  clearly  evident.  A  plot  of  the 
contours  of  constant  <t>  are  shown  in  Figure  31.  Deviations  of  the 
contours  from  straight  lines  appear  as  the  Van  der  Waals  gap  is 
approached. 

The  total  current  in  the  x-direction  Ix  is  given  by 
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(51) 


rx  -  /"  Jx 
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'  2cJo 

The  effective  resistance  R  of  the  lithium  borate  glass  sample  can 
be  written  as 


Numerical  calculations  have  been  made  of  R  as  a  function  of  the 
ratio  c/w  for  various  values  of  the  ratio  L/2w,  and  the  results 
are  plotted  in  Figure  32.  As  the  ratio  c/w  decreases  below  the 
value  of  0.1,  the  effective  resistance  increases  sharply.  This 
effect  may  be  of  some  importance  in  solid  state  batteries 
fabricated  with  thin  films. 

At  the  interface  between  the  lithium  borate  glass  and  the 
InSe,  there  can  be  a  potential  barrier  that  the  Li+  ions  must 
surmount  in  crossing  the  interface. 
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